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We describe in detail a previously published measurement of CP violation in B°-B° oscillations, 
based on an integrated luminosity of 425.7fb~ 1 collected by the BABAR experiment at the PEPII 
collider. We apply a novel technique to a sample of about 6 million B° -A D* * * § **+ £~vt. decays selected 
with partial reconstruction of the D* + meson. The charged lepton identifies the flavor of one B 
meson at its decay time, the flavor of the other B is determined by kaon tagging. 

We ^determine a CP violating asymmetry Acp = ( N(B°B° ) - N(B 0 B°))/(N(B°B 0 ) + 
N(B°B 0 )) = (0.06 ± 0.ir±g |!)% corresponding to A C p = 1 - \q/p\ = (0.29 ± 0.84±i;|f) x 10" 3 . 
This measurement is consistent and competitive with those obtained at the B factories with dilepton 
events. 

PACS numbers: 13.25.Ft, 13.20.He, 13.20.Gd 


I. INTRODUCTION 

The time evolution of neutral B mesons is governed by 
the Schrodinger equation : 

-i — V = WS> (1) 

where + ^ 2 | B°) and B° = ( bd ) and B° = 

(bd ) are flavor eigenstates. The Hamiltonian PL = M— 
is the combination of two 2x2 hermitian matrices, 
Mt = m, rt = r, expressing dispersive and absorp¬ 
tive contributions respectively. The two eigenstates of 
PL, with well-defined values of mass (rrip , mp) and decay 
width (r L , Th), are expressed in terms of B° and B°, as 

\Bl) = p\B°) + q\B°) (2) 

I B H ) = p\B°) — q\B°), 

where 


Q _ / ir^ 2 /2 

P " V M 12 -iT 12 /2' y) 

The process of B°-B° flavor mixing is therefore governed 
by two real parameters, |iU/i2 1, |P 12 1 5 and by the phase 
</>i 2 = arg(— T 12 /M 12 ). 

The value of \Mi 2 \ is related to the frequency of B°- 
B° oscillations, Am, by the relation: 

Am = m H - m L = 2\M 12 \, (4) 


*Now at: University of Tabuk, Tabuk 71491, Saudi Arabia 

TNow at: Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, F-75252 Paris, France 

tNow at: University of Huddersfield, Huddersfield HD1 3DH, UK 

§ Deceased 

^Now at: University of South Alabama, Mobile, Alabama 36688, 
USA 

**Also at: Universita di Sassari, 1-07100 Sassari, Italy 


whereas the following expression relates the decay width 
difference AT to |Ti 2 1 and <j>i 2 : 

AT = Tl — T# = 2|Fi21 cos^>i2. (5) 

A third observable probing mixing is the CP mixing 
asymmetry 

= 2(l-|2|) = Aha„ te , (6) 

where V = prob(H° -A B°) is the probability that a 
state, produced as a B°, decays as a B°, V = prob(H° —> 
B°) is the probability for the CP conjugate oscillation, 
the second equality holds if \q/p\ — 1 , and the last if 
ITia/Mial 1. 

In the Standard Model the dispersive term M\ 2 is dom¬ 
inated by box diagrams involving two top quarks. Owing 
to the large top mass, a sizeable value of Am is expected. 
The measured value Am = 0.510 ± 0.004 ps -1 (Tj is con¬ 
sistent with the SM expectation. The period corresponds 
to about eight times the B° average lifetime. 

As only the few final states common to B° and B° 
contribute to |Fi 2 |» small values of AT and Acp are ex¬ 
pected. One of the most recent theoretical calculations 
based on the SM [2f, including NLO QCD correction, 
predicts : 

A cp = 1 - \q/p\— ^cp = — (2.4±g;|) x 10~ 4 . (7) 

Sizeable deviations from zero would therefore be a clear 
indication of New Physics. A detailed review of possible 
NP contributions to CP-violation in B°- B° mixing can 
be found in [3]. In this paper, we describe the measure¬ 
ment of Acp performed by the BABAR collaboration with 
a novel technique, previously published in {¥], which, due 
to the analysis complexity, requires a more detailed de¬ 
scription. 

This article is organized as follows. An overview of the 
current experimental situation and the strategy of this 
measurement are reported in S ec. |Hl The BABAR detec¬ 
tor is described briefly in Sec. |III| Event selection and 
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sample composition are then described in Sec. IV Tag¬ 
ging the flavor of the B meson is described in Sec. [V] 
The measurement of Acp is described in Sec. m the 
fit validation is described in Sec. m the discussion of 
the systematic uncertainties follows in Sec. VIII while 
we summarize the results and draw our conclusions in 
Sec. [IX] and |Xj 


II. EXPERIMENTAL OVERVIEW AND 
DESCRIPTION OF THE MEASUREMENT 


In hadron collider experiments, bb pairs produced at 
the parton level hadronize generating the b hadrons, 
which eventually decay weakly. In B factories, pairs 
of opposite flavor B-mesons are produced through the 
process e + e~ —> T’(4S') -A BB in an entangled quantum 
state. Because of flavor mixing, decays of two B° or 
B° mesons are observed. If CP is violated in mixing, 
V ^ V and a different number of B°B° events with re¬ 
spect to B°B° is expected. The asymmetry is measured 
by selecting flavor tagged final states /, for which the 
decay B° -A f is allowed and the decay B° -A / is 
forbidden. Inclusive semileptonic decays B° —>■ £ + V(X 
have been used in the past, due to the large branch¬ 
ing fraction and high selection efficiency (unless the con¬ 
trary is explicitly stated, we always imply charge con¬ 
jugated processes; “lepton” l means either electron or 
muon). Assuming CPT symmetry for semileptonic de¬ 
cays (T{B° —► £ + V(X) = T{B° —► £~u e X)), the observed 
asymmetry is directly related to CP violation in mixing: 

M(£+£+)~ M(£-£~) _ 

Af(£+£+) + N (£-£-) AcP 

where -V(tt ) is the efficiency-corrected number of 
equal charge dilepton events after background subtrac¬ 
tion. 

Published results from CLEO [5] and the B factory 
experiments Belle E3 and BABAR 7[ [5j , based on the anal¬ 
ysis of dilepton events, are consistent with the SM expec¬ 
tation. The D0 collaboration [9], using a dimuon sample, 
obtained a more precise measurement, which however in¬ 
cludes contributions from both B° and B® mixing. They 
observe a deviation larger than three standard deviations 
from the SM expectation. Measurements based on the 
reconstruction of B° s -A D^ + £vt decays EH EE] and of 
B° -A D^ + £i7t decays [I21IT3] are compatible both with 
the SM and with D0. 

The dilepton measurements benefit from the large 
number of events that can be selected at B factories or at 
hadron colliders. However, they rely on the use of control 
samples to subtract the charge-asymmetric background 
originating from hadrons wrongly identified as leptons or 
leptons from light hadron decays, and to compute the 
charge-dependent lepton identification asymmetry that 
may produce a false signal. The systematic uncertanties 
associated with the corrections for these effects constitute 


a severe limitation of the precision of the measurements. 
Particularly obnoxious is the case when a lepton from a 
direct B semileptonic decay is combined with a lepton of 
equal charge from a charm meson produced in the decay 
of the other B. As the mixing probability is rather low, 
this background process is enhanced with respect to the 
signal, so that stringent kinematic selections need to be 
applied. Authors of m suggest that at least a part of the 
1)0 dilepton discrepancy could be due to charm decays. 

Herein we present in detail a measurement which over¬ 
comes these difficulties with a new approach. To reduce 
the background dilution from B + B~ or from light quark 
events, we reconstruct B° -A D* + £~ vg decays with a very 
efficient selection using only the charged lepton and the 
low-momentum pion (tt s ) from the D* + -A D°tt s decay. 
A state decaying as a B° (B°) meson produces £ + tt s ~ 
(£~n s + ). We use charged kaons from decays of the other 
B° to tag its flavor {K T ). Kaons are mostly produced 
in the Cabibbo-favored (CF) process B° -A DX,D -A 
K + X\ so that a state decaying as a B° (B°) meson re¬ 
sults most often in a I\ + (K ~). If mixing takes place, the 
t and the K will have the same charge. Kaons produced 
in association with the £tt s pair are used to measure the 
large instrumental asymmetry in kaon identification. 

The observed asymmetry between the number of 
positive-charge and negative-charge leptons can be ap¬ 
proximated as: 


Ag ~ Art + Acp ■ Xd, (9) 


where \d = 0.1862 ± 0.0023 E is the integrated mixing 
probability for B° mesons, and A r t is the charge asym¬ 
metry in the reconstruction of B° -A D* + £~i>\ decays. 

With the same approximations as before, the observed 
asymmetry in the rate of kaon-tagged mixed events is: 


A _ N(£+K+)-N(£~Kf) 
T N{£+K+) + N(£-Kt) 


Art + Aic + Acp, ( 10 ) 


where Ak is the charge asymmetry in kaon reconstruc¬ 
tion. A kaon with the same charge as the £ might also 
come from the CF decays of the D° meson produced with 
the lepton from the partially reconstructed side (ATr). 
The asymmetry observed for these events is: 


A _ N(£+K+)-N(£-Kn) 
R N(£+K+) + N(£-K~) 


Art + Ak + Acp ■ Xd 


( 11 ) 


Equations |9|[l0l and[TT| defining quantities computed in 
terms of the observed number of events integrated over 
time, can be inverted to extract Acp and the detector 
induced asymmetries. It is not possible to distinguish a 
Kp from a Kp in each event. They are separated on 
a statistical basis, using kinematic features and proper¬ 
time difference information. 
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III. THE BABAR DETECTOR 

A detailed description of the BABAR detector and the 
algorithms used for charged and neutral particle recon¬ 
struction and identification is provided elsewhere mm- 
A brief summary is given here. The momentum of 
charged particles is measured by the tracking system, 
which consists of a silicon vertex tracker (SVT) and a 
drift chamber (DCH) in a 1.5 T magnetic field. The po¬ 
sitions of points along the trajectories of charged particles 
measured with the SVT are used for vertex reconstruc¬ 
tion and for measuring the momentum of charged par¬ 
ticles, including those particles with low transverse mo¬ 
mentum that do not reach the DCH due to the bending 
in the magnetic field. The energy loss in the SVT is used 
to discriminate low-momentum pions from electrons. 

Higher-energy electrons are identified from the ratio 
of the energy of their associated shower in the elec- 
tromagetic calorimeter (EMC) to their momentum, the 
transverse profile of the shower, the energy loss in the 
DCH, and the information from the Cherenkov detector 
(DIRC). The electron identification efficiency is 93%, and 
the misidentification rate for pions and kaons is less than 
1 %. 

Muons are identified on the basis of the energy de¬ 
posited in the EMC and the penetration in the instru¬ 
mented flux return (IFR) of the superconducting coil, 
which contains resistive plate chambers and limited 
stramer tubes interspersed with iron. Muon candidates 
compatible with the kaon hypothesis in the DIRC are 
rejected. The muon identification efficiency is about 
80%, and the misidentification rate for pions and kaons 
is ~ 3%. 

We select kaons from charged particles with momenta 
larger than 0.2 GeV/c using a standard algorithm which 
combines DIRC information with the measurements of 
the energy losses in the SVT and DCH. True kaons 
are identified with ~ 85% efficiency and a ~ 3% pion 
misidentification rate. 


IV. EVENT SELECTION 

The data sample used in this analysis consists of 468 
million BB pairs, corresponding to an integrated lumi¬ 
nosity of 425.7fb _1 , collected at the T(iS) resonance 
(on-resonance) and 45fb _1 collected 40MeV below the 
resonance (off-resonance) by the BABAR detector 4~7j . 
The off-resonance events are used to describe the non- 
BB (continuum) background. A simulated sample of BB 
events with integrated luminosity equivalent to approx¬ 
imately three times the size of the data sample, based 
on EvtGen [T8] and GEANT4 jTS] with full detector re¬ 
sponse and event reconstruction, is used to test the anal¬ 
ysis procedure. 

We preselect a sample of hadronic events with at least 
four charged particles. To reduce continuum background 
we require the ratio of the 2 nd to the 0 th order Fox- 


Wolfram variables m be less than 0.6. We then select a 
sample of partially reconstructed B mesons in the chan¬ 
nel B° -A D*+Xt V£, by retaining events containing a 
charged lepton (£ = e, p) and a low-momentum pion (soft 
pion, 7T+) from the decay D* + -A D° 7r+. The lepton mo¬ 
mentum must be in the range 1.4 < p/- < 2.3 GeV/c 
and the soft pion candidate must satisfy 60 < p^+ < 
190MeV/c. Throughout this paper the momentum, en¬ 
ergy and direction of all particles are determined in the 
e + e _ rest frame. The two tracks must be consistent with 
originating from a common vertex, constrained to the 
beam-spot in the plane transverse to the e + e - collision 
axis. Finally, we combine pi- , p + and the probability 
from the vertex fit into a likelihood ratio variable (p ). A 
cut on ?y is optimized to reject background from other 
BB events. If more than one combination is found in an 
event, we keep the one with the largest value of p. 

The squared missing mass is: 

Admiss 2 = (-Eq, eam — Ed- — E() 2 — (p D , + pe) 2 , (12) 

where we neglect the momentum of the B° (pn 
ss 340 MeV/c) and identify the B° energy with the beam 
energy -E^eam i n the e + e _ center-of-mass frame; Ef and 
Pe are the energy and momentum of the lepton and p£>* 
is the estimated momentum of the D* . As a consequence 
of the limited phase space available in the D* + decay, 
the soft pion is emitted in a direction close to that of 
the D* + and a strong correlation holds between the en¬ 
ergy of the two particles in the B° center of mass frame. 
The D* + four-momentum can, therefore, be estimated 
by approximating its direction as that of the soft pion, 
and parameterizing its momentum as a linear function 
of the soft-pion momentum using simulated events. We 
select pairs of tracks with opposite charge for the sig¬ 
nal (^^"77) and we use same-charge pairs (t ± 'K s ± ) for 
background studies. 

Several processes where D* + and t~ originate from 
the same R-meson produce a peak near zero in the 
A4 m i ss 2 distribution. The signal consists of (a) B° -A 
D* + £~Di decays (primary); (b) B° -A D* + (mr)£~ve 
( D **), and (c) R° -A D* + t~v t , r“ —» £~v^v T . The 
main source of peaking background is due to cliarged- 
B decays to resonant or non resonant charm excitations, 
B + -A D* + (wk)£~ or to r leptons, and B -A D* + h~X , 
where the hadron (h = n, K, D) is erroneously iden¬ 
tified as, or decays to, a charged lepton. We also in¬ 
clude radiative events, where photons with energy above 
1 MeV are emitted by any charged particle, as described 
in the simulation by PHOTOS [2T]. We define the sig¬ 
nal region Ad m i ss 2 > —2 GeV 2 /c 4 , and the sideband 
-10 < Ad m i ss 2 < -4 GeV 2 /c 4 . 

Continuum events and random combinations of a low- 
momentum pion and an opposite-charge lepton from 
combinatorial BB events contribute to the non-peaking 
background. We determine the number of signal events 
in the sample with a minimum-^ 2 fit to the Ad m i ss 2 dis¬ 
tribution in the interval —10 < Ad m iss 2 < 2.5 GeV 2 /c 4 . 
In the fit, the continuum contribution is obtained from 


off-peak events, normalized by the on-peak to off-peak 
luminosity ratio, the other contributions are taken from 
the simulation. The number of events from combinato¬ 
rial BB background, primary decays and D** ((a) and 
(b) categories described previously) are allowed to vary 
in the fit, while the other peaking contributions are fixed 
to the simulation expectations (few percent). The num¬ 
ber of B° mesons in the sample is then obtained assuming 
that 2/3 of the fitted number of D** events are produced 
by B + decays, as suggested by simple isospin considera¬ 
tions. We find a total of (5.945±0.007)-10 6 signal events, 
where the uncertainty is only statistical. In the full range 
signal events account for about 30% of the sample and 
continuum background for about 15%. The result of the 
fit is shown in Fig. |T| 
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Kr from a CF decay has the same charge as the l. Equal- 
charge combinations are also observed from Cabibbo- 
supressed (CS) Kr production in unmixed events, and 
from CF Kr production. Unmixed CF Kt, mixed CS 
Kt, and CS Kr, result in opposite-charge combinations. 
Other charged particles wrongly identified as kaons con¬ 
tribute both to equal and opposite charge events with 
comparable rates. 

We distinguish Kt from Kr using proper-time differ¬ 
ence information. We define A Z = Z lec — Z tSLg , where 
Z lec is the projection along the beam direction of the Br 
decay point, and Z ta , g is the projection along the same 
direction of the intersection of the K track trajectory 
with the beam-spot. In the boost approximation |22j we 
measure the proper-time difference between the two B 
meson decays using the relation At = AZ/((3"fc), where 
the parameters /3 ,7 express the Lorentz boost from the 
laboratory to the Y(4S) rest frame. We reject events if 
the uncertainty a (At) exceeds 3 ps. 

Due to the short lifetime and small boost of the D° me¬ 
son, small values of At are expected for the Kr. Much 
larger values are instead expected for CF mixed Kt, due 
to the long period of the B° oscillation. Fig. [ 2 ] shows the 
At distributions for Kt and Kr events, as obtained from 
the simulation. To improve the separation between Kt 
and Kr, we also exploit kinematics. In the rest frame of 
the B°, the t and the D* + are emitted at large angles. 
Therefore the angle Oik between the t and the Kr has 
values close to 7 r, and cos 6^k close to -1. The correspond¬ 
ing distribution for Kt is instead uniform, as shown in 
Fig- [3] 

In about 20% of the cases, our events contain more 
than one kaon: most often we find both a Kt and a 
Kr candidate. As these two carry different information, 
we accept multiple candidate events. Using several sim¬ 
ulated pseudo-experiments, we assess the effect of this 
choice on the statistical uncertainty. 


FIG. 1: (color online): Afmiss 2 distribution for the selected 
events. The data are represented by the points with error 
bars. The fitted contributions from B° —> D* + £~ue plus 
B° —> D* + T~vt, peaking background, D ** events (1/3 from 
B° and 2/3 from B + decays), BB combinatorial, and rescaled 
off-peak events are shown (see text for details). 


V. KAON TAG 

We indicate with Kr (Kt) kaons produced from the 
decay of the D° from the partially reconstructed B° 
(Br), or in any step of the decay of the other B (Br). 
We exploit the relation between the charge of the lepton 
and that of the Kr to define an event as “mixed” or “un¬ 
mixed” . When an oscillation takes place, and the two B° 
mesons in the event have the same flavor at decay time, a 


VI. EXTRACTION OF A CP 

The measurement proceeds in two stages. 

First we measure the sample composition of the eight 
tagged samples grouped by lepton kind, lepton charge 
and K charge, with the fit to A4 m i ss 2 previously de¬ 
scribed. We also fit the four inclusive lepton samples 
to determine the charge asymmetries at the reconstruc¬ 
tion stage (see Eq. [9]). At this point of the analysis we 
use the total number of collected events. 

The results of the first stage are used in the second 
stage, where we fit simultaneously the cos Oik and At 
distributions in the eight tagged samples. 

The At distributions for BB, BB and BB events 
are parameterized as the convolutions of the theoreti¬ 
cal distributions ^(At'l©) with the resolution function 
1Z(At, At'): Gi(At) = /+“ Ti(At'\e)n(At,At')d(At'), 
where At' is the actual difference between the times of 
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decay of the two mesons and 0 is the vector of the physi¬ 
cal parameters. The decays of the B + mesons are param¬ 
eterized by an exponential function, P B + = r + e“l r+At I, 
where the B + decay width is the inverse of the lifetime 


T^ 1 = r + = 1.641 ± 0.008 ps [Ti . According to Ref. [25] , 
the decays of the B° mesons are described by the follow¬ 
ing expressions: 


J~ B o B o(At ) 

KB°B°(At ) 


B B OBo(At') 


£(A t') 


£(A t') 

£(At') 

£(A t') 

£(A t') 
To 


2(1+ r' 2 ) 



p-ColAt'l 
° 5 


cosh(APAt , /2) + 
cosh(APAt , /2) + 
cosh(APAt , /2) — 
cosh(APAt , /2) — 




cos(AmAt') — 
cos(AmAt') + 
cos(AmAt') — 
cos(AmAt') + 


Q 

P 

P 

Q 

P 

Q 

Q 

P 


(b + c) sin(AmAt') 
(b — c) sin(AmAt') 
(b — c) sin(Am,At ; ) 
(b + c) sin(AmAt / ) 


q 

P 

P 

q 


(13) 

(14) 

(15) 

(16) 


where the first index refers to the flavor of the B b and 
the second index to that of the Bt■ In Eqs. |13|fl6l 
T 0 = r B o _1 is the average width of the two B' mass 
eigenstates, A?n and AT are respectively their mass and 
width differences, r' is a parameter resulting from the in¬ 
terference of CF and doubly Cabibbo-suppressed (DCS) 
decays on the B p side, b and c two parameters expressing 
the CP violation arising from that interference. In the 
Standard Model the value of r' is rather small, 0(1%), 
b = 2 r' sin(2/3 + 7) cos<5', and c = —2 r' cos(2/3 + 7) sin<5', 
where (3 and 7 are angles of the Unitary Triangle 123] , 
and S' is a strong phase. Besides \q/p\, also Am, t b 0, 
sin(2/3 + 7), b , and c are determined as effective param¬ 
eters to reduce the systematic uncertainty. The value of 
AT is fixed to zero, and then varied within its allowed 
range [Tj when computing the systematic uncertainty. 
Neglecting the tiny contribution from DCS decays, the 
main contribution to the asymmetry is time independent 
and due to the normalization factors in Eqs. [Ts] and [T6} 

When the K T comes from the decay of the B° meson 
to a CP-eigenstate (as, for instance B° —> !)(*)/)(*)), a 
different expression applies: 

Pcp(At') = e _r °l At '[I ± S'sin(AmAt') (17) 
± Ccos(AmAf')], 

where the sign plus is used if the B b decays as a B a and 
the sign minus otherwise. This sample contains several 
components and is strongly biased by the selection cuts, 
therefore we take the values of S and C, and the frac¬ 
tion of these events in each sample (about 1%) from the 
simulation. 

The resolution function P{At, At') accounts for the 
uncertainties in the measurement of At, for the effect 
of the boost approximation, and for the displacement of 


the Kt production point from the B b decay position 
due to the motion of the charm meson. It consists of 
the superposition of several Gaussian functions convolved 
with exponentials. 

We determine A cp with two different inputs for 
Qx R {At), the At distribution for K B events, and take 
the mean value of the two determinations as the nom¬ 
inal result. As first input, we use the distribution ob¬ 
tained from a high purity selection of K B events on data, 
g^(At) HP . As second input, we use the distribution 


for K b events as predicted by the simulation. 


corrected using Eq. 18 


r, mc 

y i<R 


(At), 


GkT( At) = £®At) x 


'g%^(At)' 


g^(At) 


(18) 


HP 


To select the high purity K B sample, we require the 
lepton and the kaon to have the same charge. As dis¬ 
cussed above, this sample consists of about 75% genuine 
K Bl where the residual number of events with a Kt is 
mostly due to mixing. Therefore we select events with a 
second high-momentum lepton, with charge opposite to 
that of the first lepton. According to the simulation, this 
raises the K B purity in the sample to about 87%. Fi¬ 
nally, we use topological variables, correlating the kaon 
momentum-vector to those of the two leptons, to raise 
the K b purity in the sample to about 95%. 

Due to the large number of events, the fit complex¬ 
ity, and the high number of floated parameters, the time 
needed for an unbinned fit to reach convergence is too 
large, therefore we perform a binned maximum likelihood 
fit. Events belonging to each of the eight categories are 
grouped into 100 At bins, 25 a (At) bins, 4 cos O^k bins, 
and 5 A4 m j ss 2 bins. We further split the data into five 
bins of K momentum, px, to account for the dependen- 
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cies of several parameters, describing the At resolution 
function, the cos(fW) distributions, the fractions of Kr 
events, etc., observed in the simulation. 


Accounting for events with wrong flavor assignment 
and Kr events, the peaking B° contributions to the equal 
and opposite charge samples in each bin j are: 


G?+k+ O') = (1 + Are)(l + Ak){(1 - fZm-^)0 B OBo(j)+u;-g B o B om+f^(l-u; l+ )g KR (j)(l+x d Ac P )} (19) 
Ge-K-ti) = (! - Arl){ 1 - Ak){{ 1 - fK~)[{ 1 - u~)g S o B o(j) + u + G b ob°(j)] + / k “(1 - ^'~)Gk r U){ 1 - XdAcp)} 
Gf+K-ti) = (1 + Art){l — Ak){(1 — /if B )[(l — W )G b Obo(j) + U + Gb°B°(j)} + Ikh^^GkhU)^ + XdAcp)} 

Ge-K+ti) — (1 — Ari){l + Ak){{1 — /k r )[(1 — uj+ )GbO B oU) + w G B 0 B o{j)] + O') (1 — XdAcp)} 


where the probability density functions (PDFs) 
G b° b° ( 7 ^$o£$o(Af) and ^^o^o(At) 

are the convolutions of the theoretical distributions in 
Eqs. I3|I6 with the resolution function. The reconstruc¬ 
tion asymmetries A r e are determined separately for the e 
and p samples. Wrong-flavor assignments are described 
by the probabilities uA for Br and w ,=t for Br. They 
are different because Kr come from a mixture of D 
mesons, while Kr are produced by D° decays only. The 
parameters Ix r {pk) describe the fractions of Kr tags 
in each sample as a function of the kaon momentum. 
Due to the different charge asymmetry of the Kr and 
the Kr events, (see equations 10 and 111, the fitted 


values of f K (pk) and \q/p\ are strongly correlated. 
The J^(pk) fractions can be factorized as: 

&Q/p\) = &Q/p\ = 1) X g^(\q/p\) (20) 

where the fxiW/pl = 1) parameters are left free in the 
fit and g ±:t (\q/p\) are analytical functions. In order to 
limit the number of free parameters in the fit, the frac¬ 
tions of Kr events in the B + sample are computed from 
the corresponding fractions in the B° samples: 


fK„{B+) = f±±{\q/p\ = l)xR 


±± 


( 21 ) 


where R ±± are correction factors obtained from the sim¬ 
ulation. 

The combinatorial background consists of B + and B° 
decays with comparable contributions. A non-negligible 
fraction of B° combinatorial events is obtained when the 
lepton in B —> D*Xlv decay is combined with a soft pion 
from the decay of a tag-side D* + . As the two particles 
must have opposite charges, the fraction of mixed events 
in the B° combinatorial background is larger than for 
peaking events. In the simulation we find that the ef¬ 
fective mixing rate of the combinatorial events depends 
linearly on the kaon momentum according to the relation: 


Xd 


= Xo° mb ( a + b • Pk), 


where 


comb 
\ — 


^comb 


2(1 


( 22 ) 


(23) 


and x C omb = Am comb T^° mb . In this expression, Am comb 
and T^° mb are the mass difference and lifetime measured 
in combinatorial events. To account for this effect, we use 
for B° combinatorial background the same expressions as 
for the signal (see Eq. 19), with the replacements: 


,^comb _ r> 

y B°B° ~ y B°B° 


/->comb _ r> 

y B°B° ~ y B°B° 


y comb 

X-d 

v comb ' 
AO 


y comb 

xCd 

v comb 5 
AO 


(24) 


-»comb 
^ B°B° 


= G 


B°B° 


1 - x c J 


comb 


!-Xo 


comb 5 


1 _ v c( 

/^comb (-> _ x A d 

y B°B° ~ y B°B° 


comb 


1 ^ Xo 


comb ' 


22 


A?n' 


comb 


and 


_comb 

r B° 


^comb/ 


The parameters a and b in Eq. 
are determined in the fit. 

The probabilities to assign a wrong flavor to Br in the 
combinatorial sample are found to be different in mixed 
and unmixed events. 

Different sets of parameters are used for peaking and 
for combinatorial events, including lifetimes, frequencies 
of B° oscillation, detector related asymmetries, whereas 
the same value of \q/p\ is used. For B + combinatorial 
events, the same PDFs as for peaking B + background 
are employed, with different sets of parameters. 

The distribution t/ con t (At) of continuum events is rep¬ 
resented by a decaying exponential, convolved with a 
resolution function similar to that used for B events. 
The effective lifetime and resolution parameters are de¬ 
termined by fitting simultaneously the off-peak data. 

We rely on the simulation to parameterize the cos Oik 
distributions. The individual cos Oex shapes for the eight 
BB tagged samples are obtained from the histograms of 
the corresponding simulated distributions, separately for 
Kr and Kr events, whereas we interpolate off-peak data 
to describe the continuum. 
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The normalized At distributions for each tagged sam- butions from peaking, BB combinatorial, and continuum 
pie are then expressed as the sum of the predicted contri- background events: 


T eK (At, a A t, Mmiss 2 , cosde tK , Pk\t b o , Am, \q/p\) = 

(1 fB+ (Almiss ) /CP (Almiss ) /comb (/W miss ) /cont (Almiss 2 ))<5frc (At, a At, cos Oi,k,Pk) 

+fB+(Mmi as 2 )G?K (At, CTAt, COS 6i,k,Pk) + fcp(Mmi ss 2 )G?K (At, <7At,COs6t,K,PK) + /comb (At miss 2 ) GpK comb ( At, (JAt, COS Bi,k,Pk) 
+fcomb(Aiini sa 2 )Ge i K Comb (At, a At, cos 9(,K,Pk) + f CO r,t(Mmi a 2 )Gl°K t (At, (J At, COS0e,K ,Pk) (25) 


where the fractions of peaking B + (/ B +), CP eigenstates 
( fcp), combinatorial BB (/ com b), and continuum (/ con t) 
events in each Ad m i ss 2 interval are taken from the re¬ 
sults of the first stage. The fraction of B° (/? omb ) and 
of B+ events (/+ mb = /comb - / c ° omb ) in the combinato¬ 
rial background have been determined from a simulation. 


The functions Gfx(j) for peaking B° events are defined 
in Eq. [l9l the functions Gfx (j), Gex(j), £«f C ° mb (/)> 


G?K Conlh ti) an£ i S'I'k* (j) are the corresponding PDFs for 
the other samples. 

For a sample of B° signal events tagged by a kaon from 
the Bt meson decay, the expected fraction p of mixed 
events in each px bin depends on Am, tb o, and a; , and 
reads 


pexp _ 


o 


J i+K+ 




C T 

y e-K~ 


~,B U 

*t+K+ 


Sf-K- 


4 e+K- 


+ gfj 


( 26 ) 


i~K+ 


19 


taking 


where the functions G^x are obtained from Eq. 
into account only the contributions from B^ events. 

We estimate this fraction by multiplying the likelihood 
by the binomial factor 


dU 

/^i D T _ 

w-'ra — 


m 


N'NJ 


( pexp)AT m( i _ pexpj 




( 27 ) 


where N m and N u are the number of mixed and unmixed 
events, respectively, in a given px bin for each subsample. 
These are obtained as the sums of the numbers of mixed 
events tagged by a kaon of a given charge, 

= Nm t X+ + N mj x- ( 28 ) 

N u = N Ut x+ + N UyK - ■ ( 29 ) 

Finally, N = N m + N u . 

The corresponding value of P^ v comb for the sample 
of B° combinatorial events tagged by a kaon from the 


Bt meson decay depends on Am comb , Tg° mb , and the 
wrong flavor assignment probability for the mixed and 
unmixed subsamples. 

For a sample of B° events tagged by a kaon from the 
Bt meson decay, the expected fraction of mixed and un¬ 
mixed events, tagged by a positive charge kaon, depends 
on Acp and the detector charge asymmetries. For the 
B° signal sample this fraction reads 


pexp 

r m{u) ,K+ 


g 


e+(e-)K+ 


g 


B° 

£+(£~)K+ 


+ G 


£-(t+)K~ 


( 30 ) 


We estimate this quantity by multiplying the likeli¬ 
hood by the binomial factor 


. -yB^ _ AT m ( u )! 

m(u),K+ - N IjV . „ I 

(^(V + ) W+ 

(1 -P^ u) k + ) N ^ )k 


( 31 ) 


For a sample of B° events tagged by a kaon from the 
Bx meson decay, the fraction of mixed events depends on 
a/ . The fraction of mixed and unmixed events, tagged 
by a positive charge kaon, depends on the detector charge 
asymmetries and on Acp- Analogously, the correspond¬ 
ing fractions for a sample of B + events tagged by a kaon 
from the Bt or Br meson decay, give information on the 
detector charge asymmetries. 

The same values of Acp and Ax are shared between 
signal and combinatorial B° samples. The values of P^ xp 
and P)) 5 ) 1 / K+ for all the subsamples are obtained from 
the ratio of integrals of the corresponding observed PDFs. 
We maximize the likelihood 
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n K i x nn c ™wc+( fc ) c u( fc ) 


n =1 
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where the indices i, j, m and n denote the mixed (un¬ 
mixed) events, tagged by a kaon of a given charge, the 
index k denotes the px bin, and the index l denotes the 
signal (combinatorial background) subsample, according 
to the B meson charge (B° or B + ), and the tagging kaon 
category (K T or I< R ). 

A total of 168 parameters are determined in the fit. 
To reach the convergence of the fit, we use a three step 
approach. In the first step we fit only the parameters 
describing the B R event fractions, whereas all the other 
parameters are fixed to the values obtained on simulated 
events. In the second step we fix the B R fractions to 
the values obtained in the first step and we float only 
the parameters describing the resolution function. In the 
last step we fix the resolution parameters to the values 
obtained in the second step and we float again all the 
other parameters together with A cp- 



FIG. 2: (color online): At distributions for Kt (a) and for 
Kr (b), as predicted by the simulation. 
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FIG. 3: cos (Ock) distributions for Kt (a) and for Kr (b), as 
predicted by the simulation. 



VII. FIT VALIDATION 


Several tests are performed to validate the result. We 
analyze simulated events with the same procedure we 
use for data, first considering only B° signal and adding 
step by step all the other samples. At each stage, the 
fit reproduces the generated values of A cp (zero), and 
of the other most significant parameters (A r e, Ak , Am, 
and Tgo). 

We then repeat the test, randomly rejecting B° or B° 
events in order to produce samples of simulated events 
with A cp = ±0.005, ±0.01, ±0.025. Also in this case the 
generated values are well reproduced by the fit. By re¬ 
moving events we also vary artificially A r e or Ak , testing 
values in the range of ±10%. In each case, the input val¬ 
ues are correctly determined, and an unbiased value of 
A cp is always obtained. A total of 67 different simulated 
event samples are used to check for biases. 

Pseudo experiments are used to check the result and 
its statistical uncertainty. We perform 173 pseudo¬ 
experiments, each with the same number of events as 
the data. We obtain a value of the likelihood larger than 
in data in 23% of the cases. 

The distribution of the fit results for A cp, obtained us- 
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ing the MIGRAD minimizer of the MINUIT f35j physics 
analysis tool for function minimization, is described by 
a Gaussian function with a central value biased by 
—3.6 x 10~ 4 (0.4 cr) with respect to the nominal result. 
We quote this discrepancy as a systematic uncertainty 
related to the analysis bias. 

The pull distribution is described by a Gaussian func¬ 
tion, with a central value —0.48 ± 0.11 and RMS width 
of 1.44 ± 0.08. The statistical uncertainty, is, there¬ 
fore, somewhat underestimated. As a cross check, by 
fitting the negative log likelihood profile near the min¬ 
imum with a parabola, we obtain an estimate of the 
statistical uncertainty from the A cp values for which 
— log£ = — log£ m in + 0.5 pQ. This result is in good 
agreement with the RMS width of the distribution of the 
pseudo-experiments results, which we take as the statis¬ 
tical uncertainty of the measurement. 


VIII. SYSTEMATIC UNCERTAINTIES AND 
CONSISTENCY CHECKS 


We consider several sources of systematic uncertainty. 
We vary each quantity by its error, as discussed below, 
we repeat the measurement, and we consider the varia¬ 
tion of the result as the corresponding systematic uncer¬ 
tainty. We then add in quadrature all the contributions 
to determine the overall systematic uncertainty. 

Peaking Sample Composition: we vary the sample com¬ 
position in the second-stage fit by the statistical uncer¬ 
tainties obtained at the first stage; the corresponding 
variation is added in quadratuture to the systematic un¬ 
certainty. We then vary the fraction of B° to B + in the 
D** peaking sample in the range (50 ± 25)% to account 
for (large) violation of isospin symmetry. The fraction 
of the peaking contributions fixed to the simulation ex¬ 
pectations is varied by ±20%. Finally we conservatively 
vary the fraction of CP-eigenstates by ±50%. 

BB combinatorial sample composition: the fraction of 
B + and B° in the BB combinatorial background is deter¬ 
mined by the simulation. A difference between B + and 
B° is expected when mixing takes place and the lepton is 
coupled to the tag side tt s from B° -A D* + X decay. We 
then vary the fraction of B° to B + events in the combina¬ 
torial sample by ±4.5%, which corresponds to the uncer¬ 
tainty in the inclusive branching fraction B° -A D* + X. 

At resolution model: in order to reduce the time in the 
fit validation, all the parameters describing the resolution 
function, which show a weak correlation with \q/p\, are 
fixed to the values obtained using an iterative procedure. 
We perform a fit by leaving free all the parameters and 
we quote the difference between the two results as a sys¬ 
tematic uncertainty. 

Kr fraction: we vary the fraction of B + —> KrX to 
B° -A KrX by ±6.8%, which corresponds to the un¬ 
certainty on the ratio BR(D*° -A K~X)/BR{D* + —>■ 
K~X). 


TABLE I: Breakdown of the main systematic uncertainties on 
A cp- 


Source 

SAcp [10 ] 

Peaking sample composition 

+ 1.50 
-1.17 

Combinatorial sample composition 

± 0.39 

At resolution model 

±0.60 

Kr fraction 

±0.11 

Kr At distribution 

±0.65 

Fit bias 

+0.58 

-0.46 

CP eigenstate description 

- 

Physical parameters 

+0 

-0.28 

Total 

+ 1.88 
-1.61 


K r At distribution: we use half the difference between 
the results obtained using the two different strategies to 
describe the Kr At distribution as a systematic uncer¬ 
tainty. 

Fit bias: we consider two contributions: the statisti¬ 
cal uncertainty on the validation test using the detailed 
simulation, and the difference between the nominal result 
and the central result determined from the ensemble of 


parameterized simulations, described in Sec. VII 


CP eigenstates description: we vary the S and C pa¬ 
rameters describing the CP eigenstates by their statisti¬ 
cal uncertainty as obtained from simulation. 

Physical parameters: we repeat the fit setting the value 
of AT to 0.02 ps -1 instead of zero. The lifetime of the 
B° and B + mesons and Am are floated in the fit. Alter¬ 
natively, we check the effect of fixing each parameter in 
turn to the world average. 

By adding in quadrature all the contributions de¬ 
scribed above, and summarized in Table III we determine 


an overall systematic uncertainty of + 'g? x 10 


1-3 


IX. RESULTS 

We perform a blind analysis: the value of A cp is kept 
masked until the study of the systematic uncertainties is 
completed and all the consistency checks are succesfully 
accomplished; the values of all the other fit parameters 
are not masked. 

After unblinding we find: A cp = (0.29 ±0.84) x 10 -3 . 
We report in Table [TT] the fit results for the most signifi¬ 
cant parameters. The value of Am is consistent with the 
world average, while the value of tro is slightly larger 
than expected, an effect also observed in the simulation. 
By fixing its value to the world average, the A cp re¬ 
sult decreases by 0.18 x 10 -3 . This effect is taken into 
account in the systematic uncertainty computation. Fig¬ 
ures [3] and [5] show the fit projections for At and cos Oik, 
respectively. 

A sizeable charge asymmetry is observed at the recon¬ 
struction stage, for both e and p reconstruction and at 
the K tagging stage, somewhat smaller than that ob¬ 
served in the simulation. As the size of A r e is the same 
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TABLE II: Fit results for the most significant parameters with their statistical uncertainty. Second column: fit to the data; 
third: fit to simulated events; last: values of the parameters in the simulation at generation stage. The detector asymmetries 
in the last column are obtained from the comparison of the reconstruction efficiencies for positive and negative particles in the 
simulation. 


Parameter 

Data fit 

Simulation fit 

From MC information 

A cp 

Arec,e 

Arec,R 

Ak 

Tgo (PS) 

Am (ps -1 ) 

Arec,e (comb) 
Arec,fi (comb) 
r B o (comb) (ps) 
Am(comb) (ps -1 ) 

0.0003 ± 0.0008 
0.0030 ± 0.0004 
0.0031 ± 0.0005 
0.0137 ± 0.0003 
1.5535 ±0.0019 
0.5085 ± 0.0009 
0.0009 ± 0.0004 
0.0024 ± 0.0005 
1.3132 ±0.0017 
0.4412 ± 0.0008 

0.0003 ± 0.0005 
0.0097 ±0.0002 
0.0084 ± 0.0003 
0.0147 ±0.0001 
1.5668 ±0.0012 
0.4826 ± 0.0006 
0.0085 ± 0.0002 
0.0103 ±0.0002 
1.2898 ±0.0012 
0.4000 ± 0.0005 

0 

0.0090 ±0.0003 
0.0091 ±0.0003 
0.0151 ±0.0001 
1.540 

0.489 

0.0095 ± 0.0002 
0.0102 ±0.0002 


for the e and the /i samples, it is reasonable to suppose 
that the main source of charge asymmetry at the B° re¬ 
construction stage is due to the n s . 

Recently the BABAR collaboration published a mea¬ 
surement of the asymmetry Acp between same-sign in¬ 
clusive dilepton samples £ + £ + and i~ i~ using the com¬ 
plete recorded data set [8]. The systematic errors of the 
two analyses are essentially uncorrelated. The correla¬ 
tion of the statistical errors is estimated to be on a level 
below 10 percent. 


X. CONCLUSIONS 

We present a new precise measurement of the parame¬ 
ter governing CP violation in B°- B° oscillations. With 
a technique based on partial B° -A D* + £~i7( reconstruc¬ 
tion and K tagging we find 

Acp = 1 - \q/p\ = (0.29 ± 0.84^;g) x 10 -3 , 

where the first uncertainty is statistical and the second 
is systematic. The corresponding asymmetry, 

A CP =* 2(1 - |^|) = (0.06 ± 0.171°;!)%, 

is consistent with and competitive with the results from 
dilepton measurements at the B factories, LHCb m and 
-D0 [9]. We observe no deviation from the SM expecta¬ 
tion PJ. 


x 10‘ 

„ 4000 

a 

to 

rO 3000 
O 

'w 2000 

Q> 

'C 

1 'two 

s. x io : 

Jo 2000 

o 

\ 

in 

C 1000 

C 

0) 

0 

0.1 

0.05 

0 

-0.05 

- 0-1 


:(o) 

i I* IS* 

-□ Tog Side ft ■ K 

Peoking 

■■Tog Side • i] 

0KC j * 

’ ■ Reco Side .ft ' 

Peeking - B 

; ■ Reco 
. SKG 

:(b) 

i. . 

-(c) 

A l*K* 

1 i 

-(d) 

A l*K- 


(e) Asymmetry 

UM+ i +v. 


5 ,10 . 

. At (ps) 


,10 X 

At(ps ) 


P 




-10 -7.5 


-2.5 


2.5 


7.5 


' 0 , X 

At(ps) 


FIG. 4: (color online): distribution of At for the continuum- 
subtracted data (points with error bars) and fitted contri¬ 
butions from peaking AY, background AY, peaking Kr and 
background Kr, for (a) £ + K + events, (b) £~ K~ events, (c) 
£~ K + events, (d) £ + K~ events, (e) raw asymmetry between 
£ + K + and £~K~ events. 
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FIG. 5: (color online): distribution of cos 6i_k for the 
continuum-subtracted data (points with error bars) and fitted 
contributions from peaking Kt, background Kt, peaking Kr 
and background Kr, for (a) £ + K + events, (b) £~ K~ events, 
(c) £~ I\ + events, (d) £ + K~ events. 
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